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' (57) Abstract: A method and apparatus are provided for performing a contrast-based dynamic range management (C-DRM) al- 

j gorithm. The apparatus comprises a C-DRM processor that performs the C-DRM algorithm of the invention in order to compress 

. an image input to the C-DRM processor down to a desired gray scale range for observation on a display. The C-DRM processor 

J decouples adjustment of image mean values (low frequency) and image contrast values (high frequency), and manages mean and 
contrast separately. The use of separate mean and contrast modification functions improve on other known compression techniques 

) by providing a more deterministic behavior and reduced complexity, allowing, for example, independent management of negative 

► and positive contrasts. The C-DRM processor can also automatically adapt to the dynamic range of an input image so that the input 

* image thereby applying the minimal compression needed to display the image. 
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A METHOD AND APPARATUS FOR PERFORMING 
A CONTRAST BASED DYNAMIC RANGE 
MANAGEMENT ALGORITHM 

BACKGROUND OF THE INVENTION 

The present invention relates to a method, a memory medium and an 
apparatus implementing and performing a contrast-based dynamic range management 
("C-DRM") algorithm, for compressing an intensity dynamic range of an input image 
to a reduced intensity dynamic range supported by an available display device and, 
5 moreover, for doing so in a fashion which maximizes the displayed image contrast 
and detail. The compression is performed, in accordance with the invention, by 
directly and separately managing the mean (low frequency) and contrast (high 
frequency) content of the input image. This affords more deterministic behavior and 
reduced complexity, while enabling automatic adaptation to the dynamic range of the 

10 image, and results in minimizing the extent of compression required to enable 
displaying the image with optimized contrast on an available display device. The 
invention more particularly relates to such a method, a memory medium and an 
apparatus applicable to cardiac x-ray imaging and accordingly operable at video rates 
(typically 30 frames per second) and having low latency (i.e., a time duration from 

15 image acquisition to display of around 150 ms and thus effectively producing a real 
time display) to enable eye-hand coordination and which does not require interactive 
tuning of performance by a physician, permitting the physician to concentrate on 
performing an on-going medical procedure in reliance upon the displayed images, 
e.g., guiding a catheter through a blood vessel of a patient, without being distracted by 

20 any requirement for interactive tuning of the image acquisition, processing and 
display apparatus. 

Digital X-ray imaging is a well-known, non-contact technique for 
observing, in real time, interior aspects of an object. In practice, an X-ray beam is 
generated and targeted on the object of interest. A detecting device is positioned on 
25 the other side of the object and detects the X-rays transmitted through the object. The 
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detected X-ray signals are converted to digital signals that represent various features 
in the object, are further processed, and the resulting signal is displayed on a display 
device such as a CRT. 

One of the fundamental image processing problems in digital X-ray 
imaging is the need to transform the intensity dynamic range of the input image to the 
dynamic range supported by an available display device. Typically, the intensity 
dynamic range of an image exceeds the dynamic range of the display mechanism by 
several times. The objective of the transformation accordingly is to compress the 
"DC" or mean component of the different regions comprising the image so that the 
dynamic range (typically '256 gray levels) of an available display device may be 
utilized in a fashion which maximizes the displayed image contrast - taking into 
account, as well, the generally recognized limitation of the human eye of discerning 
only 256 gray levels. 

Common approaches to achieving such transformations are known as 
"unsharp masking" and "extended dynamic range (EDR)", the latter a special 
adaptation of the former. The conventional approach to performing the EDR 
algorithm, in general, is simply to subtract a portion of the mean from the input 
signal. However, in some situations, this approach results in important contrast 
(higher frequency) information either being removed from the image or being 
artificially enhanced and thereby introducing artifacts. In certain improved 
implementations of the EDR algorithm, however, compression is less likely to result 
in a loss of higher frequency contrast information or an introduction of artifacts. 

A more recent type of EDR function, employed in a cardiac feasibility 
study relating to a predecessor system relative to the present invention, is computed 
and implemented in the following manner: 



y(ij) = GAMMA[x(i,j)- BOOSl[x(ijj\ ] , 
Y{i,j) = GAMMA[x{i,j)-BOOST\x[ij)\ 



(1) 
(1) 
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where: 

y(i, j) is the (i,j) lh pixel value of the output image; 
x(i, j) is the (i,j) th pixel value of the input image; and 

x (i, j), i.e., (x_bar(i,j)), is the local spatial mean intensity value of 
5 the (i,j) ,h pixel, derived from a BOXCAR average. (In practice, intensity is directly 
related to X-ray count, but the relationship is rather complex. The x-rays are 
converted to photons which, in turn, are converted the /electrons, in an x-ray imager. 
The Electrons are then digitized by an analog-to-digital converter ("A/D" or "ADC") 
and basic image corrections are performed, such as adjustments to gain, offset and 
10 scaling, after which the image is ready to be processed for display.) 

A graphical representation of the EDR processing is shown in Fig. 1. 
The intensity value x of an input pixel (ij) of an image is first processed by a 
BOXCAR (moving average) function 12 that determines the local mean intensity 
value at that (ij) pixel location. (An "x" is used herein to designate an input intensity 
15 at a pixel location and, thus constitutes an individual scaler value; by contrast, an "X" 
designates the intensity image value, and thus is a vector value.) The BOXCAR 
function 12 utilizes a neighborhood of pixels, which includes and is centered on the 
input pixel, to calculate the local spatial mean intensity value x(ij) - (see, terms of 
Equation (1), supra). 

20 As illustrated in Fig. 1, BOOST LUT 14 comprises a look-up table 

(LUT) which specifies the intensity reduction of the input image signal x(i j) as a 
function of the local spatial mean intensity value *(ij). An adder (ADD) 18 
combines the (negative) output of BOOST LUT 14 and the (positive) output of LUT 
16 (see Equation (1)) and supplies the result to GAMMA LUT 20, which then 

25 compresses the result of the unsharp masking, or subtraction, operation of an ADD 1 8 
to 256 levels (8 bits per pixel, or 8 bpp). The LUTs 14 and 20 are indexed by the 
appropriate pixel intensity values given in equation (1). Thus, each of the BOOST 
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LUT 14 and the GAMMA LUT 20 jointly manages both mean and contrast 
modification functions. Thus, this more current EDR processing algorithm, while an 
improvement over previous compression/transformation algorithms (e.g., which 
merely subtracted the local mean intensity signal from the input signal), is relatively 
5 complex, yet does not permit simultaneous, independent control of the mean and 
contrast modification functions. 

Another problem with the current EDR processing algorithm is that of 
inconsistent contrast management resulting in exaggeration of negative contrast 
regions in the image. When a region in the image, such as a vessel filled with dye, for 
10 example, has an image intensity which is less than the surrounding local mean 
intensity value, EDR processing may exaggerate the negative contrast associated with 
the darker region when it subtracts the local mean intensity values from the intensity 
values associated with the darker region. This exaggerated negative contrast may 
result in artifacts, which can lead to misdiagnosis. 

15 Where an image includes multiple areas having potentially differing 

mean levels of gray, or when images of objects embedded in such areas have 
respective, different gray levels, or if respective gray levels of an object and its 
background have similar values, the contrast parameters of the display window must 
be adjusted to enhance the visibility of these differences in order to obtain a diagnosis 

20 of the underlying structure being imaged. Thus, when a viewer's attention is shifted 
from one object to another, where the contrast of one combination of object and 
background differs significantly from the contrast of another combination of object 
and background, various display window parameters relating to contrast adjustment 
must be changed. Without such adjustments, the image will appear either excessively 

25 faint or excessively bright, such that all detail critical to an effective diagnosis is 
absent. As a result, in order to obtain a diagnosis, it is often necessary, during the 
course of shifting attention among areas of differing contrast, for the physician to 
make numerous contrast adjustments to the display window. This can be disruptive 
in, e.g., radiology and mammography procedures and is altogether unacceptable in 

30 cardiac procedures, underway at the time. 
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Contrast of a cardiac digital X-ray image must be managed in a 
deterministic and consistent manner to achieve optimum results. For example, 
cardiologists perform diagnoses by examining the apparent thickness of a coronary 
artery, as revealed in X-ray imaging by a contrast medium injection. Because of 
5 X-ray physics, the artery thickness may lead to a modulation of the underlying 
background gray-level, i.e., contrast. Contrast consistency, particularly in dye filled 
vessels, is important, since processing-induced contrast changes in vessels may be 
interpreted as coronary disease leading to misdiagnosis. Consequently, it is desirable 
that artery contrast (i.e., not intensity) pass through DRM processing with 

10 deterministic and linear gain. However, minimal and deterministic modification to 
contrast can be tolerated, particularly in X-ray imaging of the lungs, where there is 
significant compression of the mean and minimal clinically relevant information. In 
fluoroscopic mode, the cardiologist is focused on the placement of interventional 
devices, and the rendering of the corresponding contrast is more directed by the 

15 visibility of these tools. In prior systems, however, contrast functions were managed 
in a non-linear manner, which varied with the local mean ~ leading to artifacts and 
confusion and resulting, in at least some cases, in increased patient exposure in efforts 
to position the X-ray tube and an image intensifier such that a vasculature of interest, 
e.g., over a spine diagram, could be viewed satisfactorily. 

20 Accordingly, it is desirable to provide a compression algorithm that 

overcomes the deficiencies of the aforementioned EDR processing algorithms and 
which permits managing the mean (low frequency) and contrast (high frequency) 
content of an image, separately and directly. Further, it is desirable to provide a 
compression algorithm that overcomes problems associated with inconsistent contrast 

25 management. Yet other limitations exist in even the more current EDR processing 
algorithm. For example, an X-ray imaging system is generally provided with 
controllable settings that allow the user to manually select one of three dynamic 
ranges. The EDR processor then subtracts a particular percentage of the local mean 
intensity value from the input image, intensity value, based on the setting selected by 

30 the user. If the user fails to select the appropriate setting which best accommodates 



5 



WO 01/69532 



PCT/US01/04918 



the dynamic range of a given image, the displayed image may have poor image 
quality. In many cases, this may result in the loss of more high frequency contrast 
information than is necessary to perform the compression. 

It accordingly is desirable to provide a compression algorithm which 
5 adaptively adjusts to the dynamic range of the image, so that high frequency contrast 
information is preserved, while applying minimal compression to display the image in 
a more deterministic manner and with reduced complexity. 

BRIEF SUMMARY OF THE INVENTION 

The invention provides a method, a related memory medium and an 
apparatus for implementing and performing a contrast-based dynamic range 

10 management (C-DRM) algorithm, to compress input image data of a broad dynamic 
range of intensities to a reduced dynamic range (e.g., 256 levels) of an available 
display device, by managing the image's mean (low frequency) and contrast (high- 
frequency) content, separately and directly. A BOXCAR (moving average) device 
takes an intensity image, as input, and produces an output image in which the 

15 intensity of each pixel represents the local mean intensity of the input image. In a 
first processing path, a local mean estimate of the intensity of a pixel of the input 
image is processed by a mean modification function to produce a modified local mean 
intensity of the pixel. In a second path, independent of the first path, the intensity of 
the corresponding pixel is processed in relation to the local mean estimate thereof to 

20 produce a contrast estimate, the local mean estimate is processed with a contrast 
modification function to produce a result. The result then is combined with the 
contrast estimate to produce a modified contrast estimate which then is combined 
with the modified local mean to produce a modified output intensity for the 
corresponding pixel in the reduced intensity dynamic range. While a BOXCAR 

25 moving average is used in the disclosed embodiments herein, the specific method 
used for computing the local mean estimate is not essential to the C-DRM processing. 
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In a further embodiment, C-DRM is implemented so as to be adaptive 
and thereby to follow variations in image intensities while panning, but without 
causing flicker or delay modifications in the resulting display. In yet further 
embodiments, function modification of contrast and of mean is implemented by 
5 respective multiple look-up tables (LUTs), the specific table for each being selected, 
based on a maximum intensity value for the image as a whole. However, the use of 
multiple LUTs for adaptive range compression is not essential to the central concept 
herein of separately and directly managing mean and contrast components of the 
image data. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 Fig. 1 is a flow diagram of an EDR algorithm processor of the prior 

art. 

Fig. 2 is a schematic flow diagram of an X-ray imaging system 
implementing and performing a C-DRM algorithm in accordance with the invention. 

Fig. 3 is a flow diagram of an apparatus implementing and performing 
15 the C-DRM algorithm, in accordance with a first embodiment of the invention, for 
processing X-ray image signals and transforming a dynamic range thereof for display 
on an available display unit having a lesser dynamic display range. 

Fig. 4 is a flow diagram of an apparatus implementing and performing 
the C-DRM algorithm in accordance with a second embodiment of the invention. 

20 Fig. 5 is a flow diagram of an apparatus implementing and performing 

the C-DRM algorithm in accordance with a third embodiment of the invention. 

Fig. 6 is a flow diagram of an apparatus implementing and performing 
the C-DRM algorithm in accordance with a fourth embodiment of the invention. 

Figs. 7 and 8 are flow diagrams of alternative fifth and sixth 
25 embodiments of apparatus implementing and performing the C-DRM algorithm. 
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Fig. 9 is a block diagram of an alternative, parallel processing system, 
relative to Fig. 2, implementing and performing a C-DRM algorithm in accordance 
with the invention and including any of the first through sixth embodiments of Figs. 3 
through 8, respectively. 

5 Fig. 10 is a flow diagram of a seventh embodiment of the invention, 

comprising an apparatus automatically adapting the C-DRM algorithm to a variable 
input dynamic range of the input image data, which is usable in any of the foregoing 
embodiments of Figs. 3 through 8 and may be realized in the hardware architecture 
presented in Fig. 9. 

DETAILED DESCRIPTION OF THE INVENTION 

10 Fig. 2 is a flow diagram of an x-ray imaging system 1 which may 

implement and perform the C-DRM processing algorithm of the invention. The 
system 1 comprises a table 4 on which an object 6 which is to be imaged (e.g., a 
human body) rests so as to be irradiated by X-rays emitted by an X-ray tube 8, driven 
by a generator 9 under exposure management control by a control and preprocessing 

15 computer 10. X-rays passing through the object 6 are attenuated in differing amounts 
depending upon the character of the body portions, and conditions therein, of the 
object 6 being imaged. Thus, the variously attenuated X-rays impinge upon and are 
detected by a digital detector 5, the image signal output of which is supplied to the 
control and preprocessing computer 10, all in conventional fashion. As is likewise 

20 conventional, the computer 10 controls various functions of the X-ray imaging 
process including, e.g., offset, gain, bad pixel correction, exposure management (e.g., 
as to the generator 9 and X-ray tube 8) and scaling operations and, further, 
preprocesses the received X-ray image data so as to produce a corrected X-ray image 
data output X. 

25 A display optimizing computer 12 includes a memory 14 which stores 

the image data, X, output by computer 10, which is to be processed in accordance 
with the C-DRM algorithm of the invention by the computer 12 and, further, performs 

8 



WO 01/69532 



PCT/US01/04918 



related display optimizing processing, such as spatial and temporal filtering, zoom, 
and display control functions, for producing an optimized output signal Y. Signal Y 
is supplied to, and displayed by, display unit 16. The imaging system 1 of Fig. 2 is 
adapted to perform C-DRM algorithm processing in accordance with the present 
5 invention and including each of the embodiments of Figs. 3-8 and 10, as hereinafter 
discussed. As used herein, "adapted to" and the like refer to processing components 
(e.g., programmable digital processor, application specific integrated circuit (ASIC) 
and the like) programmed or arranged to process data to provide a desired output. 

A feasibility system, as shown in Fig. 2, was implemented using 
10 Datacube™ MAXPCI ASIC's whose operations were specified in software. The 
software was Imageflow™, a Datacube™ language, used to setup the ASIC's and 
related control and data flow paths. 

As seen in Fig. 1 (prior art) and as described above, conventional 
algorithms, collectively called "EDR", use GAMMA and BOOST in a "composite 

15 function" form: GAMMA(x(ij) - BOOST(x_bar(ij))). In other words, the value at 
which GAMMA is evaluated actually depends, itself, on the value of BOOST at 
x_bar. This composite function couples the effects of GAMMA and BOOST and 
makes the algorithm less "tractable," i.e., it is harder to isolate out, or extract, the 
individual effects of changes to GAMMA and BOOST. 

20 The invention overcomes this problem by circumventing the composite 

function formulation that is a defining feature of the EDR algorithms. The two 
modification functions F_c and F_m are not composed with each other— i.e., the 
argument of each is only M_x, which comes from the BOXCAR. The values of F_c 
and F_m can be combined in the lookup tables ("LUT's") of each of the embodiments 

25 of Figs. 3 to 7 without defeating the decoupling, because neither F_c nor F_m is 
evaluated at an argument that includes the other function. 

As a result of this decoupling, each data attribute (i.e., as to each of 
contrast and mean) may be modified directly, instead of through a chain of nested 
functions. Each modification function also concerns only one attribute, either 

30 contrast or mean. The modification of the attributes is therefore performed separately 
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of, instead of jointly with, the other. Note that in EDR, the only way to adjust 
contrast at x_i,j, for example, is to change both GAMMA and BOOST jointly— 
because y_i j is defined by the composite 
(GAMMA()...+BOOST0) 

Fig. 3 illustrates a block diagram of a C-DRM processor of the display 
optimizing computer 12 of Fig. 2, in accordance with a first embodiment of the 
invention. Each block of the C-DRM processor represents one or more functions, 
where: 

X = input intensity 

M x = local mean estimate 

C x = contrast estimate 

F M () = mean modification function of M x 

F c () = contrast modification function M x , and with respect to 
embodiments of Figures 6 through 8, of M x , 

Y = C-DRM output intensity 
M Y = modified local mean 
C Y = modified contrast 

The C-DRM processor computes for Y as follows: 

Y = M Y (1+C Y ), where: (2) 

M Y =F M (M X ), (3) 
C Y =F c (Mx)*C x ,and (4) 
C X =(X-M X )/M X (5) 
As is apparent from equations (3) and (4), the modified local mean 
(My) and the modified contrast (C Y ) are separately and directly managed, 
respectively, in producing the C-DRM output intensity Y, in equation (2). F_m, the 
mean modification function, may achieve a linear or non-linear transformation of the 
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mean component of the image.. F_c, the contrast modification function, is restricted 
to applying a linear gain to the input contrast to construct the output contrast. 

In accordance with the implementation of the first embodiment of the 
invention shown in Fig. 3, the input intensity value X of an input pixel in the image 
is first processed by the "BOXCAR" (or moving average) function 30 to produce the 
local mean estimate M x . The specific method for computing the moving average is 
not essential to the C-DRM method. The BOXCAR function 30 utilizes a 
neighborhood of pixels, typically next-adjacent to and including the input pixel being 
processed, to calculate the local mean estimate The BOXCAR function 30 acts 
as a low pass filter for separating structures of interest from variations in the 
detection output which are less relevant, e.g., gradually varying background. In 
experimental implementations, a 32x32 BOXCAR resulted in ringing (that is, an 
intensity oscillation) along contrast field vessels contrast filled vessels, heart to lung 
transition, heart to diaphragm, contrast filled catheter, etc. A 64x64 BOXCAR avoids 
that deficiency by spanning a sufficient spatial area, or extent, such that contrast- 
-liquid filled coronaries do not unduly influence the mean estimate of the 
background and sudden background changes, such as heart to lung, are more 
gradual; accordingly, a 64x64 BOXCAR is preferred. 

The local mean M x is then inputted to arithmetic unit 31, which 
processes X in accordance with the equation (5) above, and outputs the contrast 
estimate C x to multiplier 33. Contrast LUT 32 processes M x in accordance with the 
contrast modification function F c ( ) of equation (4) above and supplies the contrast 
gain, the output F C (M X ) thereof, to multiplier 33 which performs the multiplication 
function of equation (4) and the produces the modified contrast C Y output. 

The mean LUT function 34 processes the local mean estimate M x by 
the mean modification function F M ( ) of equation (3), above, to output the modified 
local mean M Y . The mean modification function F M () compresses the low frequency 
component of the input intensity X, directly and separately from (i.e., independently 
of) the contrast component. Adder 36 then performs the function equation (2), above, 
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outputting the C-DRM output intensity Y to the display 11, thereby to reproduce the 
image thereon in the reduced intensity dynamic range supported by the display 11. 

The following specifications were employed in an actual reduction to 
practice of the invention: 

5 MEAN FUNCTION MODIFICATION 

DRM shall compress the background of the input image (in the range 0 
to NxABD) into an 8 bits component (in the range of 0 to 255) by applying a 
monotonically increasing modification function My=F M (Mx) depending upon Mx 
only, defined by: 

10 F M (0)=0 

F M (ABD)=100 

F M (NxABD)<255 

The monotonically increasing function F M means that F M (X+l)> F M 

(X)VX. 

15 As an example of the mean modification function, Fm() transforms the 

mean from being linearly related to X-ray counts to being linearly related to an 
equivalent water thickness represented by the X-ray counts. This is achieved through 
the application of well known mathematical models of X-ray physics. The 
transformation from "intensity" to "thickness" is commonly constrained to map the 

20 assumed mean heart intensity, i.e., the automatic brightness detection "ABD" value, 
to a desired output mean, typically 100. The assumed input mean heart intensity is 
achieved through the conventional exposure management and rescaling functions of 
the X-ray system 1 of Fig. 2, and particularly the computer 10. 

CONTRAST FUNCTION MODIFICATION 



25 



It is assumed here that the input image is in a linear encoded format. 
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The contrast modification function C Y = F c (Mx)*C x is implemented, 
preferably, so as to comply with the requirement that F c shall be monotonically 
increasing, which means that: 

F c (C+l , Mx) > F c (C,Mx) VC, Mx. (5 A) 

5 The contrast modification function Cy is a linear gain which is a 

function of M x , as seen in the above equations (and accordingly is subject to variation 
along with variation in the mean). In practice, the gain is nominally unity in the heart 
and in other, relatively darker anatomical regions where clinically relevant 
information exists. The contrast gain must smoothly transition to values less than 
10 unity in the relatively bright lung region, to suppress the detail, which is not clinically 
relevant and may constitute a possible source of distraction, and to minimize clipping 
as the mean component of the resulting image approaches the maximum intensity of 
255. In other words, the available contrast afforded by the maximum intensity of 255 
must be carefully allocated to clinically relevant data. 

15 The C-DRM data flow diagrams in Fig 3 may be optimized for 

real-time processing in accordance with the further, alternative embodiments of Figs. 
4 and 5. 

Fig. 4 illustrates a second embodiment of the invention, partially 
optimized for real-time processing relative to the embodiment of Fig. 3. Specifically, 

20 the Fig. 3 division function 31 and the multiply function 35 are performed within the 
LUT component 42 of Figure 4, thereby reducing the Fig. 3 computation functions by 
one division function and one multiply function. In Fig. 4, input intensity X is 
processed by BOXCAR function 40, as in BOXCAR 30 of Fig. 3, and the resulting 
local mean estimate M x is then directed to adder 41 for subtractive combination with 

25 X, and which produces X-M x as an output. M % is also supplied to 1/M Contrast Mean 
LUT function 42, the output of which is supplied to multiplier 43 and the product 
output of which, C Y *M Y , representing a combined, modified contrast and mean result, 
is then directed to a second adder 45. M s furthermore is supplied to the mean LUT 
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function 44 which outputs, as does mean LUT 36 in Fig. 3, the modified local mean 
M Y . M Y is combined in adder 45 with the output of multiplier 43, above described, to 
produce the compressed C-DRM output intensity signal Y=M Y (1+C Y ). 

A third embodiment of the invention, which is fully optimized for real- 
-time implementation, is illustrated in Fig. 5. The computation of the invention is 
further simplified by eliminating the ADD 41 and unsigned (that is, without 
arithmetic sign) inputs to the MULTIPLY 43 of Fig. 4. BOXCAR 50, as with the 
BOXCARS 12, 30 and 40 of the prior embodiments, performs a moving average on 
the input intensity X to output the local mean estimate M x , Mx is supplied to LUT 51, 
where it subjected to the therein-designated modification function, and the resulting 
product then is directed to multiplier 52 for multiplication with X, the resulting 
product being supplied to adder 54. M x is also supplied to LUT 53, where it is 
processed by the modification function therein indicated i.e., the mean modification 
function is multiplied by the negative complement of the contrast modification 
function (1-F 0 (M X )) and the resultant output is supplied to adder 54, for additive 
combination with the output of multiplier 52, and adder 54 then outputs the 
compressed, C-DRM output intensity Y. 

Due to the movement of computational elements of prior embodiments 
(e.g., Figs. 3 and 4) into LUTs 51 and 53, the relationship of LUTs 51 and 53 to the 
contrast and mean modification functions is complex. Nevertheless, the processing 
functions of Fig. 5 will be readily understood from the following, which leads to the 
definition of the C-DRM output intensity Y: 

Y=CyMy+My, where: (6) 

Cy=fcCx and 

fc(Mx)=fc 

Y=fcCxMy+My, where: (7) 
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_ x- Mx 
CX= Mx 



Y=fc(X-Mx)(My/Mx)+My (8) 

r-^u-m^ (9, 



My=fm(Mx)=fm 



CIO 

Y - X ^- Fm ^ (13) 



r-*(4£) + >(l-J&) (14) 



As can be seen from the succession of equations (6) through (14), the 
function modifications for contrast and mean are each handled directly and separately 
of the other, even though concurrently processed in two respective and different 
primary paths. All functions with the exception of multiplication 52 and addition 54 
in this embodiment of Fig. 5 are implemented in the LUTs 51 and 53, optimizing the 
speed of processing and thus affording a real-time implementation. 

The embodiments of the invention shown in Figs. 3-5, when applied to 
digital cardiac images, were experimentally found, in common, to provide a clinically 
accurate but aesthetically suboptimal rendering of the lung field. Input positive 
contrasts (x > Mx) were at risk of being clipped to 255 unless a very small contrast 
15 
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gain was applied but, as a result thereof, "contouring" or "cartooning" resulted, 
which is also undesirable. Note that for negative contrasts (x<Mx), x is less than Mx, 
thereby bounding Cx between 0-100%. For positive contrasts (x>Mx), x may be 
many times greater than Mx; therefore, Cx is more loosely upper bounded. This 
5 positive contrast issue provides challenges in the lung portion of an image, where 
there is significant local variation in contrast. This challenge is addressed in 
accordance with the invention by managing positive and negative contrasts 
separately. In this way, negatively contrasted detail, such as contrast-fluid (e.g., 
iodine) filled vessels and inserted clinical devices (e.g., catheters) in bright areas, 
10 such as the lungs, may pass through C-DRM with strong gain and be readily visible to 
the clinician. Conversely, positive contrasts in brights areas preferably are 
suppressed, thereby avoiding aesthetically undesirable rendition of bright areas, such 
as lung fields. 

As background relevant to the embodiment of Fig. 6, in cardiac X-ray 
15 imaging, the bright lung field is generally unimportant clinically, but nevertheless it is 
highly desired that it have a natural appearance. The contrast and mean of the lungs 
may be strongly compressed, ensuring the majority of the display dynamic range is 
allocated to anatomical areas of clinical interest but while maintaining a natural 
appearance. However, contrast-liquid filled coronaries and clinical devices may be 
20 imaged over the lungs as a result of the projection of certain imaging angulations. In 
this situation, the contrast of the clinically relevant information must be preserved and 
not compressed. To achieve the above objectives and in accordance with Fig. 6, the 
sign Sex of the contrast estimate Cx (Scx=X-Mx) is measured and Sex is used to 
manage positive and negative contrasts, separately achieving the dual objectives 
25 stated above. 

More particularly, Fig. 6 depicts the fourth embodiment of the 
invention in which positive and negative contrasts are managed separately, by 
supplying the sign, S c „ of the input contrast, or contrast estimate, as an input to the 
contrast modification LUT 62. The physical size of contrast LUT 62 is doubled, so as 
30 to accommodate the additional (sign) bit for the input address. The remainder of Fig. 
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6 corresponds to Fig. 3. The BOXCAR 60 performs a local mean estimate in the 
same manner as BOXCARS 12, 30, 40, and 50. With X and Mx as inputs, function 
66 calculates the input contrast (i.e., contrast estimate) Cx in the same manner as 
function 31 of Fig. 3. Sex, the sign of Cx, and Mx are inputs to the contrast 
5 modification function Fc(Scx, Mx) LUT 62. Mx is modified to My via the mean 
modification LUT 64. My is then input to the MULTIPLY 65 and to the ADD 69. 
The ADD 69 produces the final output Y. 

The embodiment of Fig. 6 may be computationally simplified for 
realization in the same manner as the first embodiment of Fig. 3 was simplified, in the 
10 successive second and third embodiments of Figs. 4 and 5, respectively. The only 
difference is the need to provide Sex as an input to the contrast modification function 
Fc() of LUT 62 Equations (6)-(14) similarly apply. 

The fifth and sixth embodiments of the invention, Figs. 7 and 8 
respectively, comprise two alternative, computationally simplified embodiments of 
15 Fig. 6 and which parallel the relationship of Figs. 4 and 5 to Fig. 3. Figs. 7 and 8 
tradeoff LUT size versus the need for signed multiplication components, or functions, 
and are detailed to the level required for implementation on dedicated hardware, such 
as an ASIC or a digital signal processor DSP. 

The relationship of the fifth embodiment of the invention, of Figure 7, 
20 to Fig 6 parallels the relationship of Fig. 5 to Fig. 3. 

Fig. 7 is labeled with bit size designations relevant to practical 
implementations of apparatus implementing and performing the C-DRM algorithm of 
the Fig. 7 flow diagram. Thus, image input X is designated "unsigned 12 bpp", i.e., 
an unsigned 12 bit per pixel input image signal. Similarly, bit depths of intermediate 
25 processing results are likewise labeled in bits per pixel ("bpp"). The BOXCAR 70 
performs the same operation as performed by BOXCARS 12, 30, 40, 50, and 60 of 
the earlier embodiments. Component 72 subtracts Mx from X to determine the 1 bpp 
sign of Cx, Sex. Sex is concatenated with Mx in a concatenator operator (i.e., a bit 
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extractor/combiner) 73 to form a composite (unsigned) 12 bpp address Mx, by which 
to address LUTs 74 and 76. 

Since Mx represents the mean of the low frequency component of the 
image data, it is reasonable to drop the least significant bit of Mx from the LUT 
5 address in order to reduce LUT size, which reduces memory requirements. The 
present embodiment illustrates only a reduction by 1 bit, which reduction has no 
perceivable impact on the computation fidelity of the results. However, 
experimentation has shown that several least significant bits ("lsb's") of Mx may be 
dropped from the LUT address before quantization errors become readily visible; 
10 hence, the illustration of a signal bit reduction is not limiting. 

The respective computations performed in LUTs 74 and 76 are the 
same as those performed in LUTs 51 and 53, respectively, of Fig. 5, with Fc now 
addressed by Sex and Mx. 

LUT 74 computes: 

15 l/Mx*Fc(Scx, Mx)*Fm(Mx) (15) 
LUT 76 computes: 

Fm(Mx)*( 1 - Fc(Scx, Mx)) ( 1 6) 

As with multiplier 52 in Fig. 5, multiplier 78 of Fig. 7 outputs the 
multiplication product of X and the output of LUT 74. The output of LUT 74 has a 

20 fractional component and is implicitly fixed point. The BIT EXTRACTOR 77 
extracts the integer portion of the output result of multiplier 75 and discards the 
fractional component. The ADD 78 performs an addition of the outputs of LUT 76 
and bit extractor 77 and produces the final result Y. LUT 79 clips the output of LUT 
78 to ensure an 8 bit result (i.e., computationally, the output of LUT 78 may be 

25 greater than 255 and thus is clipped by LUT 79 to ensure the 8 bit result). 
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Fig 8, the sixth embodiment of the invention, represents an alternative 
realization to that of Fig. 7, relative to the fourth embodiment of Fig 6. The 
relationship of Fig. 8 to Fig. 6 parallels that of Fig. 4 to Fig. 3. Fig. 8 requires less 
memory for LUTs than Fig. 7 due to the reduced size of LUT 86 relative to LUT 76. 
5 This reduction in LUT size comes at the cost of the additional signed multiply 
function 88 in Fig. 8. Elements 84-89 of Fig. 8 correspond to elements 74-79 of Fig. 
7, respectively, subject to the above differentiation between these embodiments in 
that Fig. 8 uses less cache but is slightly more computationally intensive since 
processing the extra sign bit, as will be appreciated by a comparison of the "signed" 
10 designations in Fig. 8 compared to the "unsigned" designations in Fig. 7. Fig. 8, 
moreover, has a combiner 81 for combining the sign bit S cx with the unsigned 11 bpp 
output of cathenator 83. Further, in Fig. 8: 

LUT 84 performs: 1/M X F C (S CX ,M X )F M (M X ) (17) 
LUT 86 performs: F M (M X ) (18) 

15 The functions in accordance with Fig. 7 were performed in a 

multiprocessor system, as shown in Fig. 9 employing, by way of example and not 
limitation, Mercury hardware using four (4) Motorola PowerPC G3 processors 90, 92, 
94, and 96 in parallel as the computational device with Mercury Raceway as data 
flow communication device. Image data from the detector 15 in a matrix of 

20 1024x1024 pixels, 12bpp, is progressively scanned in four (4) equal amounts, 
indicated in Figure 9 by Direct Memory Addressing (DMS) of lines 1-256 to memory 
9, DMA lines 257-512 into memory 93, DMA lines 513-778 into memory 95, and 
DMA lines 777-1024 into memory 97. This addressing format enables processing of 
data in parallel by the respective processors 90, 92, 94 and 96 in accordance with the 

25 C-DRM algorithm, as hereinabove explained, (i.e., Figs. 3 to 8). Processor assembler 
software and Mercury MCOS software is utilized to instruct the hardware devices to 
perform data communications in coordination with appropriate control signals and 
CDRM computation on the data in coordination with data availability. The BOXCAR 
function (illustrated as items 12, 30, 40, 50, 60, 70, and 80 in Figures 1 and 3-8 

30 respectively) is decomposed to two serial one dimensional averages, vertical followed 



WO 01/69532 



PCT/US01/04918 



by horizontal. Since the image is partitioned by lines among the four processors 90, 
92, 94, 96, the vertical boxcar inputs are shared among the neighboring processors on 
the vertical boundaries 98a, 98b and 98c. Due to the moving average computation of 
the vertical boxcar, the result must be propagated from 90 to 92 to 94, and finally to 
96 with the final result passed back up, in each instance, as shown by the 
double-headed arrows on the vertical boundaries 98a, 98b and 98c. As an alternative 
embodiment to the use of the aforesaid vertical boundaries, a corresponding 
replication of boundary lines in the data transfer input to each of the memories 91, 93, 
95 and 97 may be employed, consistent with the moving average BOXCAR function 
of the processors. Once the CDRM computation has been completed, the image result 
remains partitioned in 91, 93, 95, and 97 and is propagated from the four memories to 
a next image processing stage or to display. 

During feasibility clinical trials of implementations, and performances, 
of the prior art EDR Fig. 1 algorithm, several observations were made regarding 
perceived needs for improvement in aspects of the use and performance which have 
been addressed in the C-DRM in accordance with the invention. As one factor, it had 
been contemplated that a technician would tune the C-DRM to the dynamic range of 
the imagery during a cardiac x-ray acquisition; this, however, was considered a 
distraction by the clinician and was therefore avoided. As a result, EDR "high" was 
always selected, which assumed a maximum intensity dynamic range of the system 
(e.g. 15*ABD value in the feasibility clinical trials). However, the true dynamic 
range of the imagery could be significantly less. When the dynamic range of the 
imagery was significantly less than 15*ABD, EDR compressed the data to a subset of 
the 255 values available. This resulted in imagery of sub-optimal image quality in 
terms of brightness, contrast, and visually perceivable bit starvation. Offline tuning 
of EDR demonstrated that the impact was readily perceivable, thus leading to an 
interest in having the C-DRM mean and contrast modification functions 
automatically adapt to the dynamic range actually present in the input imagery rather 
than assume the maximum range numerically possible. 
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As background relevant to a discussion of the embodiment of Fig. 10, 
the intensity dynamic range of cardiac X-ray images varies significantly, depending 
on field of view, gantry (tube/imager) angulation relative to the patient, and lung 
content in the image. As is common knowledge, if an objective is to compress an 
5 input dynamic range of [0, 600] to a display range of [0, 255] employing a 
compression algorithm and if, instead, one assumes erroneously an input dynamic 
range of [0, 4000], the compression algorithm will produce a suboptimical result. 
Such a hypothetical, or exaggerated, value of the input dynamic range actually can 
exist, or be interpreted as existing, in an input image by virtue of anomalies or 
10 defects, e.g., as produced by erroneous bright pixels caused by imager defects. 

To overcome those problems and, in accordance with the embodiment 
of the invention shown in Fig. 10, the C-DRM algorithm automatically adapts to the 
intensity dynamic range of the image, affording a significantly broader dynamic range 
management enabling one to take full advantage of the actual dynamic range of the 

15 imagery than typically is accomplished than when manual adjustments are made, e.g., 
by a technician, for tuning the C-DRM "strength" during a cardiac x-ray acquisition. 
Automatic adaptation likewise avoids possible distraction of the clinician; likewise, it 
avoids the undesired results, where manual tuning is not done correctly, of a reduced 
dynamic range of the imagery and sub-optimal image quality in terms of brightness, 

20 contrast and bit starvation. Such automatic adaptation, moreover, should be of low 
latency (at rates of frames, not seconds) and must be smooth, so the imagery does not 
abruptly flicker during the adaptation. This feature is referred to hereinafter as the 
"MAX" function, in that the C-DRM algorithm adapts to a temporally low pass 
filtered maximum value of M K , as a measure of the image's intensity dynamic range. 

25 The MAX feature thus measures the dynamic range of the input image and then tunes 
the compression algorithm accordingly. 

Particularly, the MAX measurement is spatially low pass filtered 
(output of BOXCAR 100) thereby to suppress erroneous bright pixels, such as imager 
defects, in the load mean estimate Mx. Mx is then temporally filtered in MAX 108 
30 (e.g., a one dimensional BOXCAR or an infinite impulse response filter) to provide 
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dynamic range continuity in the output Mx', thereby suppressing rapid changes in the 
dynamic range of the imagery and, e.g., suppressing flicker in the displayed image 
when the clinician is panning. 

As shown in Fig. 10, the BOXCAR function 100 calculates the local 
5 mean estimate M x value in the manner discussed in relation to Fig. 3. The local mean 
estimate M x is then further processed by MAX function 108 to determine a maximum 
M x of the local mean estimate M x , which is then supplied to each of the contrast LUT 
function 102 and the Mean LUT function 104 , each corresponding to a particular 
sub-range of the total dynamic range. This function basically determines the dynamic 

10 range over the whole image; that information then is output as M x ' and used in LUT 
addressing. The object of M x ' is to impose only the amount of compression necessary 
- e.g., some small amount of compression of flat imagery and substantial compression 
of dynamic imagery, thereby increasing the performance. (Without the MAX 
function, a common amount of compression is applied to all imagery, regardless of 

15 the dynamic range actually present.) M x ' accordingly selects a LUT from a set of 
LUTs, and the selected LUT is then indexed by the usual Mx. 

For example, assume that the maximum intensity dynamic range of a pixel of 
an image, prior to compression, is from 0 to 1500. Each of the lookup tables in the 
Contrast LUT 102 and in Mean LUT 1044 will correspond to a respective sub-range 

20 of the total dynamic range. In one practical device, contrast LUT 22 and mean LUT 
24 each comprises 16 lookup tables which together cover the total dynamic range 
(i.e., the total dynamic range preferably is divided into 16 sub-ranges) Mx_prime 
(M x ) is an additional input to the LUTs. However, the invention is not limited with 
respect to the number of lookup tables utilized by the Contrast LUT 22 and in Mean 

25 LUT 24 

The local mean estimate M'x is supplied to the contrast LUT function 102, 
which processes Mx and Mx' by the contrast function FC(Mx, Mx') and outputs to 
multiplier 103. Mx* is also input to the mean LUT function 104, which outputs 
Fm(Mx, Mx'), the adaptively modified local mean MY', to multiplier 105 and ADD 
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106. Multiplier 105 performs the function of equation (4) to output the adaptively 
modified product Cy'*MY', which ADD function 106 combines with MY' in 
accordance with equation (2) to produce the adaptive C-DRM output intensity Y\ 

The contrast and mean values of the look-up tables, 102 and 104 respectively 
5 of Fig. 10, are similar to those of LUTs 32 and 34 in Fig. 3, the first embodiment of 
the invention, the difference being that LUTs 102 and 104 are addressed by Mx and 
Mx' thereby further adapting their respective outputs to the dynamic range present in 
X as measured by 108, i.e., Mx'. If the dynamic range is relatively large (i.e., portions 
of the image are very bright while others are very dark), the look-up indexed in the 

10 Mean Modification LUT function 104, Fm(Mx, Mx'), results in a relatively large 
compression of the local mean intensity value and, similarly, the Contrast 
Modification LUT 102, Fc(Mx, Mx'), strongly compresses contrast in the brights to 
avoid excessive clipping. Conversely, if the dynamic range is relatively small (i.e., 
the image is relatively dark), the look-up table indexed in the Mean Modification 

15 Function 104 causes modest change to the means and the Contrast Modification LUT 
function 102 better preserves a maximum contrast in the upper portion of the dynamic 
range input signal. 

It will be understood by those skilled in the art the MAX adaptation invention 
in Fig. 10 can readily be applied to the C-DRM embodiments represented in Figs. 3-8 
20 having equivalent performance benefit at the cost of more memory to store the larger 
LUTs now additionally addressed by Mx'. 

While the invention has been disclosed hereinabove in the environment of a 
digital signal processing system, it is equally applicable to an analog system, such as 
the LC+ Cardiac X-ray product produced by the General Electric. Company and which 
25 implements the EDR algorithm discussed hereinabove. In such an analog adaptation 
and with reference to Fig. 4, the input intensity X represents (i.e., rather than a pixel, 
as in the digital system described hereinabove) the analog voltage waveform for a line 
of data, successive units X representing respective, successive lines of data separated 
by horizontal synchronization signals. In the analog implementation, BOXCAR 40 
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may comprise a simple one dimensional low pass filter. ADD 41 outputs the 
difference waveform, X-Mx, by inverting Mx and summing the waveforms. Mx is 
applied to LUT 42, and the resultant output thereof modulates the X-M x signal in 
multiplier 43 to produce an output comprising an analog waveform C y ( M y , 
5 corresponding to one line of the image. LUT 44 represent analog circuitry 
performing the mean modification function F m () on Mx, resulting in M Y . ADD 45 
sums the analog output waveforms of 43 and 44, resulting in the output analog 
waveform Y, representing an analog waveform of the reduced intensity range and 
corresponding to the input analog waveform for a line of data of the image. 

10 It will be understood by those skilled in the art that the invention has been 

described with reference to the preferred embodiments, but that the invention is not 
limited to these embodiments. Those skilled in the art will understand that 
modifications can be made to the embodiments discussed above which are within the 
scope of the invention. 
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WHAT IS CLAIMED IS: 

1 . A method of compressing an intensity dynamic range of an 
input image to a reduced intensity dynamic range of an image display signal, 
comprising: 

a) defining plural units of each frame of the input image; 

b) determining (30) a local mean estimate (Mx) of the intensity (X) of 
each of the plural units of the input image; 

c) processing each local mean estimate (Mx), separately and directly, 
in accordance with: 

i) a mean modification function (34: F M ()) to produce a 
corresponding mean modification processing result (M Y ), and 

ii) a contrast modification function to produce a 

corresponding contrast modification processing result (C y ); and 

d) combining (35, 36) the corresponding mean modification function 
and contrast modification function results (C Y ) to produce, for each local mean 
estimate (X), a respective output intensity (Y), within the reduced intensity range, for 
a corresponding unit of the image display signal. 

2. A method as recited in claim 1, wherein: 

each unit corresponds to a pixel of the input image as detected by a 
digital detector. 

3 . A method as recited in claim 1 , wherein: 

each unit corresponds to a scan line of an analog image detector. 

4. A method as recited in claim 1 , wherein the processing of each 
local mean estimate, separately and directly, further comprises: 



multiplying the contrast modification function (F c ())by the inverse of 
the local mean estimate. 
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5. A method as recited in claim 4, wherein the processing of each 
local mean estimate, separately and directly, further comprises: 

multiplying the contrast modification function (F c ())by the mean 
modification function (M Y ) and producing the multiplication product as a first 
combined, corresponding mean modification and contrast modification processing 
result; and 

the combining step further comprises combining the combined (C 
Y +M Y ), corresponding mean modification function processing result (M Y ) and contrast 
modification function processing result (M Y ) with the mean modification processing 
result (M y ) to produce, for each local mean estimate, the respective output intensity 
(Y), within the reduced intensity range, for the corresponding unit of the image 
display signal. 

- 6 . A method as recited in claim 1 , wherein the processing of each 
local mean estimate, separately and directly, further comprises: 

multiplying the mean modification function by a negative complement 
of the contrast modification function and producing the multiplication product as a 
first combined, and corresponding, mean modification function and contrast 
modification function processing result; and 

multiplying the contrast modification function by (5 1) the mean 
modification function and by an inverse of the local mean estimate (X) and producing 
the multiplication product as a second combined, and corresponding, mean 
modification function and contrast modification function processing result; and 

the combining step further comprises combining the first and second 
combined, and corresponding, mean modification function and contrast modification 
function processing results to produce, for each local estimate, a respective output 
intensity (Y), within the reduced intensity range, for the corresponding unit of the 
image display signal. 
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7. A method as recited in claim 1, further comprising: 

determining (66) a contrast estimate (Cx)as a function of a ratio, of a 
difference between the intensity (X) of a unit of the input image and the local mean 
estimate (Mx) thereof to the local mean estimate (M x ); 

5 determining (66) a sign of the contrast estimate (S cx ); and 

controlling the processing of each local mean estimate (Mx) by the 
contrast modification function (62) by employing the sign of the contrast estimate 
(S cx ) to manage positive and negative contrasts separately. 

8. A method as recited in claim 1 , further comprising: 

10 processing (108) the local mean estimate (My by temporally filtering 

same and producing a resultant, modified local mean estimate having enhanced 
dynamic range continuity, substantially independent of rapid changes in the dynamic 
intensity range of the input image; and 

controlling the processing of each local mean estimate (Mx), including 
15 the processing thereof by the mean (104) and contrast (102) modification functions, in 
accordance with the modified local mean estimate (Mx). 

9. A method as recited in claim 8, further comprising: 

processing each local mean estimate (Mx) and the mean (F M ()) and 
contrast modification function processing results and the output intensity (Y) 
20 corresponding to each unit of the input image, as respective, multiple bit values; and 

limiting (77, 87) the number of bits of the output intensity of the 
corresponding unit of the image display signal to a predetermined number of bits, 
consistent with display range capabilities of an available display device. 

1 0. A method as recited in claim 9, further comprising: 
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performing the processing of each local mean estimate in accordance 
with the mean and contrast modification functions and the multiplication of the 
modification results of same for producing multiplication products, in accordance 
with accessing respective look-up tables (74, 76 and 84, 86) from which are produced 
the corresponding multiplication products and wherein each of the local mean 
estimate, the mean modification function and the contrast modification function 
processing results and the respective multiplication products thereof constitute 
respective, multiple bit values. 

11. A method as recited in claim 10, further comprising: 

determining (72, 82) a sign value of the contrast estimate, as a 
respective, single bit value for each local mean estimate; and 

combining (73, 81) the sign of a contrast estimate with the local mean 
estimate for accessing each of the look-up tables. 

12. A method as recited in claim 10, further comprising: 

determining (72, 82) a sign of the contrast estimate, as a respective, 
single bit value for each local mean estimate; 

forming (70, 30) the local mean estimate as an unsigned number of 

bits; 

combining (73) the sign of the contrast estimate with the unsigned 
number of bits representing the local mean estimate for accessing a first look-up table 
(74) performing the contrast modification function; and 

accessing a second look-up table (86) performing the mean 
modification function using the unsigned number of bits representing the local mean 
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13. An apparatus compressing an intensity dynamic range of an 
input image to a reduced intensity dynamic range of an image display signal, each 
frame of the input image comprising plural units, the apparatus comprising: 

a filter (30) determining a local mean estimate (Mx) of the intensity 
(X) of each of the plural units of the input image; 

a computer (12) processing each local mean estimate (Mx), separately 
and directly, in accordance with: 

a mean modification function (34:FM()) to produce a 
corresponding mean modification processing result (M Y ), and 

a contrast modification function (32: F c ()) to produce a 

corresponding contrast modification processing result (C Y ); and 

the computer combining the corresponding mean modification 
function and contrast modification procession function results to produce, for each 
local mean estimate (X), a respective output intensity (Y), within the reduced 
intensity range, for a corresponding unit of the image display signal. 

14. An apparatus as recited in claim 13, further comprising: 

a digital image detector (12), each unit corresponding to a pixel of the 
input image as detected by the digital detector. 

15. An apparatus as recited in claim 1 3 , further comprising: 

an analog image detector (12), each unit corresponding to a scan line 
of the input image as detected by the analog image detector. 

16. An apparatus as recited in claim 13, wherein the processing of 
each local mean estimate, separately and directly, by the computer (12) further 
comprises: 
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multiplying the contrast modification function (F c ()) by the inverse of 
the local mean estimate(C Y ). 

17. An apparatus as recited in claim 16, wherein the processing of 
each local mean estimate, separately and directly, by the computer (12) further 
comprises: 

multiplying the contrast modification function by the mean 
modification function (M Y ) and producing the multiplication product as a first 
combined, corresponding mean modification and contrast modification processing 
result; and 

combining the combined (C Y +M Y ), corresponding mean modification 
function processing result (M Y ) and contrast modification function processing result 
(M Y ) with the mean modification processing result (M Y ) to produce, for each local 
mean estimate, the respective output intensity (Y), within the reduced intensity range, 
for the corresponding unit of the image display signal. 

18. An apparatus as recited in claim 1 3 , wherein the processing of 
each local mean estimate, separately and directly, by the computer (12) further 
comprises: 

multiplying (53) the mean modification function by a negative 
complement of the contrast modification function and producing the multiplication 
product as a first combined, and corresponding, mean modification function and 
contrast modification function processing result; 

multiplying the contrast modification function by (51) the mean 
modification function and by an inverse of the local mean estimate (X) and producing 
the multiplication product as a second combined, and corresponding, mean 
modification function and contrast modification function processing result; and 

combining the first and second combined, and corresponding, mean 
modification function and contrast modification function processing results to 
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produce, for each local estimate, a respective output intensity (Y), within the reduced 
intensity range, for the corresponding unit of the image display signal. 

19. An apparatus as recited in claim 1 3 , wherein the processing by 
the computer (12) further comprises: 

5 determining (66) a contrast estimate (QJas a function of a ratio, of a 

difference (X-M x ) between the intensity (X) of a unit of the input image and the local 
mean estimate (M x ) thereof to the local mean estimate (Mx); 

determining (66) a sign of the contrast estimate (S cx ); and 

controlling the processing of each local mean estimate (Mx) by the 
10 contrast modification function (62) by employing the sign of the contrast estimate 
(S cx ) to manage positive and negative contrasts separately. 

20. An apparatus as recited in claim 13, wherein the processing by 
the computer further comprises: 

processing (108) the local mean estimate (Mx) by temporally filtering 
15 same and producing a resultant, modified local mean estimate having enhanced 

dynamic range continuity, substantially independent of rapid changes in the dynamic 
intensity range of the input image; and 

controlling the processing of each local mean estimate (M^, including 
the processing thereof by the mean (104) and contrast (102) modification functions, in 
20 accordance with the modified local mean estimate (M x '). 

21. An apparatus as recited in claim 20, wherein the processing by 
the computer further comprises: 

processing each local mean estimate (M c ) and the mean (F M ()) and 
contrast (F c ()) modification function processing results and the output intensity (Y) 
25 corresponding to each unit of the input image, as respective, multiple bit values; and 
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limiting (77, 87) the number of bits of the output intensity of the 
corresponding unit of the image display signal to a predetermined number of bits, 
consistent with display range capabilities of an available display device. 

22. An apparatus as recited in claim 21 , wherein the processing by 
5 the computer further comprises: 

performing the processing of each local mean estimate in accordance 
with the mean and contrast modification functions and the multiplication of the 
modification results of same for producing multiplication products, in accordance 
with accessing respective look-up (74, 76 and 84, 86) tables from which are produced 
10 the corresponding multiplication products and wherein each of the local mean 
estimate, the mean modification function and the contrast modification function 
processing results and the respective multiplication products thereof constitute 
respective, multiple bit values. 

23 . An apparatus as recited in claim 22, where the processing by 
1 5 the computer (12) further comprises : 

determining (72, 82) a sign value of the contrast estimate, as a 
respective, single bit value for each local mean estimate; and 

combining (73, 81) the sign of a contrast estimate with the local mean 
estimate for accessing each of the look-up tables. 

20 24. An apparatus as recited in claim 24, where the processing by 

the computer further comprises: 

determining (72, 82) a sign of the contrast estimate, as a respective, 
single bit value for each local mean estimate; 

forming (70) the local mean estimate as. an unsigned number of bits; 
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combining (73) the sign of the contrast estimate with the unsigned 
number of bits representing the local mean estimate for accessing a first look-up table 
(74) performing the contrast modification function; and 

accessing a second look-up table (76) performing the mean 
5 modification function using the unsigned number of bits representing the local mean 
estimate. 

25. A memory medium controlling a computer to compress an 
intensity dynamic range of an input image to a reduced intensity dynamic range of an 
image display signal, by: 

10 a) defining plural units of each frame of the input image; 

b) determining (30) a local mean estimate of the intensity (Mx) of each 
of the plural units of the input image; 

c) processing each local mean estimate (Mx), separately and directly, 
in accordance with: 

15 i) a mean modification function (34:F M ()) to produce a 

corresponding mean modification processing result (M Y ), and 
ii) a contrast modification function (32) to produce a 
corresponding contrast modification processing result (C Y ); and 

d) combining (35, 36) the corresponding mean modification function 
20 and contrast modification function results (Cx) to produce, for each local mean 

estimate (X), a respective output intensity (Y), within the reduced intensity range, for 
a corresponding unit of the image display signal. 

26. A memory medium as recited in claim 25, wherein: 
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each unit corresponds to a pixel of the input image as detected by a 
digital detector. 

27. A memory medium as recited in claim 25, wherein: 
each unit corresponds to a scan line of an analog image detector. 

28. A memory medium as recited in claim 25, further controlling 
the computer to process each local mean estimate, separately and directly, by: 

multiplying the contrast modification function (F c ()) by the inverse of 
the local mean estimate. 

29. A memory medium as recited in claim 28, further controlling 
the computer to process each local mean estimate, separately and directly, by: 

multiplying the contrast modification function (F()) by the mean 
modification function (M Y ) an d producing the multiplication product as a first 
combined, corresponding mean modification and contrast modification processing 
result; and 

the combining step further comprises combining the combined (C Y - 
M Y ), corresponding mean modification function processing result (M Y ) and contrast 
modification function processing result with the mean modification processing result 
(M Y ) to produce, for each local mean estimate, the respective output intensity (Y), 
within the reduced intensity range, for the corresponding unit of the image display 
signal. 

30. A memory medium as recited in claim 25, further controlling 
the computer in accordance with: 

multiplying (53) the mean modification function by a negative 
complement of the contrast modification function and producing the multiplication 
product as a first combined, and corresponding, mean modification function and 
contrast modification function processing result; 
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multiplying the contrast modification function by (5 1) the mean 
modification function and by an inverse of the local mean estimate (X) and producing 
the multiplication product as a second combined, and corresponding, mean 
modification function and contrast modification function processing result; and 

5 the combining step further comprises combining the first and second 

combined, and corresponding, mean modification function and contrast modification 
function processing results to produce, for each local estimate, a respective output 
intensity (Y), within the reduced intensity range, for the corresponding unit of the 
image display signal. 

10 3 1 . A memory medium as recited in claim 25, futher controlling 

the computer in accordance with: 

determining (66) a contrast estimate (C x ) as a function of a ratio, of a 
differenc (X-M x )e between the intensity (X) of a unit of the input image and the local 
mean estimate (Mx) thereof to the local mean estimate; 

15 determining a sign of the contrast estimate (S cx ); and 

controlling the processing of each local mean estimate (Mx) by the 
contrast modification function (62) by employing the sign of the contrast estimate 
(S cx ) to manage positive and negative contrasts separately. 

32. A memory medium as recited in claim 25, further controlling 
20 the computer in accordance with: 

processing (108) the local mean estimate (Mx) by temporally filtering 
same and producing a resultant, modified local mean estimate having enhanced 
dynamic range continuity, substantially independent of rapid changes in the dynamic 
intensity range of the input image; and 
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controlling the processing of each local mean estimate (M*), including 
the processing thereof by the mean (104) and contrast (102) modification functions, in 
accordance with the modified local mean estimate (M x '). 

33. A memory medium as recited in claim 32, further controlling 
5 the computer in accordance with: 

processing each local mean estimate (M x ) and the mean (F M ()) and 
contrast (F c ()) modification function processing results and the output intensity (Y) 
corresponding to each unit of the input image, as respective, multiple bit values; and 

limiting (77, 87) the number of bits of the output intensity of the 
1 0 corresponding unit of the image display signal to a predetermined number of bits, 
consistent with display range capabilities of an available display device. 

34. A memory medium as recited in claim 32, further controlling 
the computer in accordance with: 

performing the processing of each local mean estimate in accordance 
15 with the mean and contrast modification functions and the multiplication of the 
modification results of same for producing multiplication products, in accordance 
with accessing respective look-up tables (74, 76, and 84, 86) from which are produced 
the corresponding multiplication products and wherein each of the local mean 
estimate, the mean modification function and the contrast modification function 
20 processing results and the respective multiplication products thereof constitute 
respective, multiple bit values. 

35. A memory medium as recited in claim 34, further controlling 
the computer in accordance with: 

determining (72, 82) a sign value of the contrast estimate, as a 
25 respective, single bit value for each local mean estimate; and 
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combining (73, 81) the sign of a contrast estimate with the local mean 
estimate for accessing each of the look-up tables. 

36. A memory medium as recited in claim 34, further controlling 
the computer in accordance with: 

determining (72, 82) a sign of the contrast estimate, as a respective, 
single bit value for each local mean estimate; 

forming (70) the local mean estimate (M x ) as an unsigned number of 

bits; 

combining (73) the sign of the contrast estimate with the unsigned 
number of bits representing the local mean estimate for accessing a first look-up table 
(74) performing the contrast modification function; and 

accessing a second look-up table (76) performing the mean 
modification function using the unsigned number of bits representing the local mean 
estimate. 
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